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We are going to discuss ion-ion and deuteron - nucleus RHIC data and show that 



> 

■ they support, if not more, the idea of the new QCD phase: colour glass condensate 

OO . 

O . with saturated parton density. 

o 

P^; INTRODUCTION 
D , 

^ . We used to think that nucleus-nucleus interaction is so complicated that the distance 



between the experimental data and underlying microscopic theory: QCD, is so large that 
we are not able to give any interpretation of the data based on QCD. My goal is to show 
that this wide spread opinion is just wrong. I hope to convince you that the nucleus-nucleus 
collisions provide such an information on the initial stage of the process, parton formation 
(see Fig. which could be (and even has been) very essential in the disscussion of the new 
QCD phase: colour glass condensate (CGC). 



The main prediction of the CGC_is_t 

saturated reaching a maximal value 



le fact that the parton density in CGC region is 



l|, Q, 3| • The space-time picture of the QCD saturation 
is shown in Fig. |21 Let us consider a virtual photon, with virtuality Q, in the rest frame 
(Bjorken frame) where the photon is the standing wave which interacts with the parton 
(colour dipoles) with size r 1/Q. In the beginning of our process we have a small number 
of partons of this size but at late time the number of partons steeply increases. At some 
moment of time the partons start to populate densely in the proton and filled the whole 
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proton disk (see Fig. EI). 

The estimate t of the value of a new scale: saturation scale which relates to the size of 
the parton when the partons started to populate densely (the critical curve in Fig. E)) we 
introduce the packing factor: 

P.F. = n = X p ^ ^ X (1) 

The saturation scale is the solution to the equation: 

n{Qs{x)) = 1 (2) 



THEORY STATUS 

The scope of this talk does not allow us to discuss theory of the high parton density systen 
but, I firmly believe, that a reader should know, our approach is based on the first principles 
of the microscopic theory (QCD). Despites rather complicated technique, that we have to use 
approaching this regime, the theoretical ideas, which we based upon, are very transparent 
and can be easily explained and digested. For the diluted system of partons the main process 
is the emission of gluon and this process leads to the famous DGLAP evolution equation. 
However, when the density of partons increases the processes of recombination, which are 
proportional to the square of density (p^), should enter to the game. The competition 
between emission (oc p), which increases the number of partons, and recombination (oc 
p^), which diminishes this number, results in the equilibrium density. The phenomenon 
of approaching the maximal density we call 'parton density saturation' and the phase of 
QCD with saturated density is the colour glass condensate. The evolution equation which 

FIG. 1: The space-time pic- 
ture of the ion-ion collision 
with the main stages of the 
evolution in the final state: 
parton formation , thermal- 
ization, Quark-Gluon Plasma 
(QGP) and hadron gas. 
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FIG. 2: Space-time picture for the deep inelastic scattering. 
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FIG. 3: The distribution of partons in the transverse plane. 



describes the saturation phenomenon is a non-hnear equation 

HQS 

which final form was 

found by Bahtsky and Kovchegov . This equation is so beautiful that I decided to present 
it here despites the lack of room. 
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where the sizes of dipoles (see Fig. H}. 
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FIG. 4: The Balitsky-Kovchegov equation. 

Eq. shows that the degrees of freedom at high energy are the colour dipoles rather than 
quark and gluon themselves Q]. The physical meaning of is the dipole-target amplitude. 
The equation says that the change of this amplitude from ytoy+dy where y = \nt = \n{l / x) 
and X is the Bjorken variable fir DIS. is equal to the probability for the incoming dipole 
to decay into two dipoles. These two dipoles can interact with the target separately and 
simultaneously. The simultaneous interaction should be taken with the negative sign which 
refiects the shadowing effect or accounting for the the recombination processes. 

The Balitsky-Kovchegov equation at the moment is the main theoretical tool for in many 
applications of the CGC dynamics despites it's very approximate nature. It plays a role of 
the mean field approach in this problem. 

It is very important to understand that we have a more general approach than the mean 
field one. This approach is based on the space-time structure of the high energy interaction 
in QCD (see Fig. El). The idea of this approach is very transparent. Let us start with 
emission of the parton at time shown by the first dotted line in Fig. |S| All parton with high 
energies were created a long before this moment. The main idea j3| is that these partons can 
be treated classically while the parton emitted at this moment can be described by QCD. 
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Moving the moment of time (see the second dotted hne in Fig. El we should include the gluon 
produced in (tj, Zi) into the parton system which we consider classically, but a new gluon 
emitted in (tj+i, -Zj+i) shall be treated in full QCD. Since the description should be the same 
in these two moments (tj and tj+i we have a constraint which leads to the equation. The 
realization of this program is rather complicated as well as technique that is required to 
understand the resulting equation but the physical basis for the JIMWLK equation 0] is 
very simple. 
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FIG. 5: The space-time structure of high energy interaction in QCD and the renormalisation 
Wilson group approach (JIMWLK equation) to high parton density QCD. 

The last remark that I would like to make in this brief theoretical introduction is related 
to the theoretical input which we will use in our description of the experimental data. It 
turns out that we need to know the solution to the non linear equation only in vicinity of 
the saturation scale together with the general understanding of the behaviour of the dipole 
amplitude in the saturation domain. Since the unitarity itself leads to — 1 the different 
equations can only specify the character of approaching this maximal value. The description 
of the RHIC data do not depend on the details of this approaching. The behaviour near 
to the saturation boundary can be found just knowing the solution to the linear evolution 
equation. 
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COMPARISON WITH THE DATA: GENERAL APPROACH. 

I believe that the most important thing is to stipulate clearly what kind of assumptions 
we are making applying the general theory to an interpretation of the experimental data 
at accessible energies and typical distances reached by the experiment. Indeed, the theory 
was formulated for the dense parton system while the experimental data exist rather in 
the transition domain where density is not very high but not low as well. In our KLMN 
approach 0, S H U, Q we used three main assumptions which we will have to refor- 



mulate at the end of the talk. 

1. At X — Xq ~ 10 we have the McLerran-Venugopalan model for the inclusive 
production of gluons with the saturation scale QKxq) ~ IGeV"^ j|3|; 

2. The RHIC region of x ~ 10^^ is considered as the low x region in which 
as ln(l/x) ~ 1 while as <^ 1. This is not a principal assumption, but it makes 
the calculations much simpler and more transparent; 

3. We assume that the interaction in the final state does not change significantly the 
distributions of particles resulting from the very beginning of the process. For hadron 
multiplicities, this may be a consequence of local parton hadron duality, or of the 
entropy conservation. Therefore multiplicity measurements are extremely important 
for uncovering the reaction dynamics. However, we would like to state clearly that wee 
do not claim that interaction in the final state are not important. We rather consider 
the CGC as the initial condition for the interaction in the final state. 



COMPARISON WITH THE DATA: HERA DATA. 



We would like to make three statements about deep inelastic scattering data from 



HERA: 



(i) models that incorporate the parton density saturation (see for example Refs. 



13 
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16l |) are able to describe the HERA data; (ii) the solution to the Balitsky-Kovchegov equation 
leads to the good description of the experimental data on the DIS structure functions 17. isl]: 
and (iii) the gluon structure function extracted from the fit of the experimental data is so 
large that the packing factors reaches the value of unity. 

The set of figures illustrate this our point of view. Fig. IB] and Fig. [7| show the comparison 
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with the experimental data the simple Golec-Biernat and Wuesthoff model in which the 
dipole-proton cross section is written in the form 

with _R^(x) = 1/Ql{x) , where 

Qlix) = Qli-Y' (5) 

The parameters were chosen from the fit of the experimental data and they are: cxo = 
23.03 mb; A = 0.288 ; xo = 3.04 10"^ Ql = IGeV^. 
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FIG. 6: j*p total cross sections in the Golec- 
Biernat and Wuesthoff model 
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FIG. 7: DIS diffraction in the Golec-Biernat 
and Wuesthoff mode . 



Fig. IHl gives the example how the Balitsky-Kovchegov equation is able to describe the 
deep inelastic structure function. 

Fig. ini gives the packing factor calculated in one of the parametrization for the extracted 
gluon structure function. 
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FIG. 8: F2 as result of solution to the Balitsky-Kovchegov equation. Picture is taken fron Ref. |l7l | 

FIG. 9: Packing factor k (see 
Eq. (HI)) from the HERA data. 

MULTIPLICITY AT RHIC 

Multiplicity is the most reliable test of our approach since the third assumption is correct 
due to entropy conservation. The physical picture of ion-ion collision in the CGC approach 
is given by Fig. EH In central region of rapidities we have a collision of two dense system 
of partons while in the fragmentation regions one system of partons is rather diluted while 
the second one turns out to be more dense than at 77 = 0. 
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FIG. 10: Ion-ion collision in the CGC approach. 



To characterize the ion- ion coUision we use several observables, namely, centrality cut, 
which related too typical value of the impact parameters; number of participants (Npart) , 
which shows the number of proton taking part in interaction; and the number of collisions 
(Ncoii) ■ All these observables are discussed in details in Nardi's talk at this conference and I 
will assume that you are familiar with them. Here, I only want to formulate the simple rule: 
if a physical observable oc Npart this observable can be described by soft physics while 
observables oc Ncoii are certainly related to typical hard processes, which can be treated 
in perturbative QCD. Soft physics at scale larger than 1 GeV is a great surprise. In simple 
words the saturation in CGC is a claim that soft physics starts at pt < Qs{x) and Qs could 
be large at high energies. 

The key relation in the CGC is that the saturation scale is proportional to the number 
of participants or 



In the case of the ion-ion collisions we have actually two saturation scales in colliding 
ions. In the simple Golec-Biernat and Wuesthoff model they have the following dependence 
on energy W and rapidity y: 




(6) 



and 



Sa-QI{A) oc Npart 



(7) 
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max 



W \ 

Ql{A;W = Wo;y = 0) {—J 



(9) 



In our papers we use the following formula for the inclusive production 



E- 



da A-KNr 



pt 



(10) 



where Xi^2 = {Pt/ ^/s) exp{^y) and ^9^1,712(3;, ) is the unintegrated gluon distribution of a 
nucleus ( for the case of the proton one of (pA should be replace by (pp.) This distribution is 
related to the gluon density by 



xG{x, 



dkl^{x,kl 



(11) 



The multiplicity can be calculated using the following formula: 

dN _ I f ^ 2 f^dcx\ 4:7TNcas 1 
dy S J * I d^p I 
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2 
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X2GA2{x2,Pt) XiGA^ixi,Pt] 



(12) 



S Jo pf 

where we integrated by parts and used Eq. (fTT|. In the KLMN treatment we assumed a 
simplified form of xG, namely, 



xG{x]pI] 



Sp] (1 - xf Pt < Qs{x) ■ 



(13) 



^5QK^) (1 - Pt > Qs{x) ; 
Inserting Eq. (jl3p into Eq. (|12j) we obtain a simple formulas: 
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^QCD / \ ^ \Vs ,max / V 

This simple formula with the saturation scales determined by Eq. (jH} and Eq. describes 
quite well the rapidity, energy and Npart dependence of the multiplicity (see Fig. [TTl Fig. [T21 
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21 



22 



23|. 



and Fig. IT^ that have been measured at RHIC 

Using Eq. p4|) we are able to predict the multiplicity distribution at the LHC energy 24 1 
The systematic errors in our predictions mostly stem from the uncertainties in the energy 
behaviour of the saturation scale (see Ref. j2J] for details). The predictions are shown in 
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FIG. 11: Rapidity dependence of measured multiplicity. 



Fig. and Fig. ^| One can see that we expect rather low value of the multiplicity in 
comparison with the majority of other approaches. 

The same formula could be generalized to describe the deuteron-nucleus collisions jll|. 
In Fig. Eland Fig. El we show our predictions for the multiplicity distribution for deuteron- 
gold collisions. Fig.lT^was published before RHIC data and one can see that our predictions 
do not agree with the experimental quite well. Therefore, we had a dilemma: the approach 
is not correct or some phenomenological parameters were chosen incorrectly. Fortunately, 
the second is the case. Indeed, we found that we have to change several parameters. It turns 
out that the number of participants calculated by us (see Nardi's talk) does not agree with 
the number of participants that experimentalist calculated using the Glauber Monte Carlo. 
In new comparison we use the experimental value for Npart- Second, we treated incorrectly 
the events with the number of participants in the deuteron less than 1. The third change 
was in the value of the saturation momentum of the proton. We have to take it just the 
same as in the Golec-Biernat and Wuesthoff model while in our old fit we took it by 30% 
less. After these corrections the comparison is not bad (see Fig. [T7j) but we cannot describe 
the fragmentation region of nucleus. 



Npart SCALING AT RHIC FOR ION-ION COLLISIONS 

One of the most interesting results from RHIC, I believe, is so called Npart scaling. As 
you can see in Fig. ITHl experimental data show that dN/dyd?pt oc Npart at all values of pt up 
to pt ~ 5 GeV. As we have discussed this fact is a strong indication that even at sufficiently 
large momenta the mechanism close to the 'soft physics' works. This scaling behaviour has 
a very simple explanation in the CGC approach, which is based on three observations jlol |: 
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FIG. 12: PHENIX data for Npart dependence 
of the multiplicity distributions. Our calcula- 
tions are labeled by 'KLN' in this figure. 
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FIG. 13: PHOBOS data on energy and Npart 
dependence of the multiplicity distributions. 
Our calculations are labeled by 'KLN satura- 
tion model' in this figure. A determines the 
energy and rapidity dependence of the satu- 
ration scale ( see Eq. (jSJ and Eq. ©) 



1. Geometrical scaling behaviour in the wide region outside of the saturation 
domainH:0(a;,P?) = (l){Q^,{x)/pf); 



2. The anomalous dimension of gluon density is ~ | for Ql{x) < pi < -M- 



(x) 



Npart ' 



3. For wide range of distances the saturation scale for a nuclear target oc Npart 
Al 

The behaviour of function looks as it is shown in Fig. □ Using such 0-s we obtain 
for ion-ion collisions: 
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(15) 
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FIG. 15: Comparison our prediction with 
other models available on the market (see Ref. 
25|. 



FIG. 14: The KLN prediction for the LHC 
energy. The full curves correspond the Golec- 
Biernat and Wuesthoff model for the energy 
dependence of the saturation scale. The dot- 
ted line describe the prediction if energy de- 
pendence of the saturation scale is determined 

n 

by the running QCD coupling. See Ref. [2^ 
for more details. 

More precise calculation shown in Fig. UHl lead to Npart scaling and can describe the 
experimental data. 



Npart SCALING AND SUPPRESSION IN DEUTERON-NUCLEUS COLLISIONS 

Using the behaviour of functions shown in Fig.nfor deuteron- nucleus collision we obtain: 

= 7^. - SoQ.WQM - (16) 

Eq. pfi|l means that we expect a suppression in comparison with wide spread opinion that 
dN/dud?pt should be proportional to Npart The experiment shows that we do not 

have such a suppression in the central and nuclear fragmentation regions of rapidities but 
it certainly exists for forward region (see Fig. |^ and Fig. 1^^ . The first result needs an 
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FIG. 16: Our prediction for the deuteron-gold collision at RHIC. 
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FIG. 17: Our final fit of the data for the deuteron-gold collision at RHIC. The left figure for 
BRAUMS data while the right one for the PHONOS data 

explanation the second one is a great success of the CGC approach. 

It should be stressed that the fact that the ration at = and at low is much smaller 
than 1 itself is a strong argument in favour of the CGC since in ordinary (Glauber) approach 
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FIG. 18: Npart scaling for gold-gold collision 
(PHOBOS data) 
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FIG. 19: Npart scaling in the CGC approach 
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FIG. 20: The behaviour 
of functions (f) for nucleus 
and proton in the CGC ap- 
proach taking into account 
the anomalous dimension. 
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it should be equal to unity. As far as I know there is no other explanation of this suppression. 

I think that the situation with a suppression is clearly illustrated by Fig. „ In this figure 
one can see two clearly separated regions with quite different physics. This separation sug- 
gests the way out: for the central and nuclear fragmentation region the anomalous dimension 
of the gluon density is not essential while in the forward fragmentation we see the effect of 
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FIG. 22: BRAHMS data in forward region 
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FIG. 23: Two different re- 
gions: Cronin enhancement 
due to interaction in the final 
state: a strong indication of 
the quark-gluon plasma pro- 
duction; and a quantum evo- 
lution region predicted bu the 
CGC approach Q| . 



quantum evolution. 



THREE ASSUMPTIONS: A NEW EDITION 



However, the fact that we do not have a suppression in the range of rapidities in central 
and nuclear fragmentation region, makes incorrect our beautiful explanation of the Npart 
scaling for ion-ion collisions. We have to assume that this scaling stems not from initial con- 
dition for our evolution but rather from strong interaction in the final state which suppress 
the production of the hadron from inside the nucleus. Only production from the nucleus 
surfaces can go out of the interacting system and can be observed. In this case we also have 
a Npart scaling. Why we see this scaling up to pt ~ 5 GeV ? We have not reached a clear 
understanding why. 

I hope that you understand that we have to change our three main assumptions under 
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the press of the experimental data. Their new edition looks as follows: 

1. At X > a;o ~ 10"^ we have the McLerran-Venugopalan model for the inclusive 
production of gluons with the saturation scale Q^(xo) ~ 1 GeV^; 

2. X < xq ^ 10~^ is low i.e. a5ln(l/a;) ~ 1 while < 1; 

3. The CGC theory determines the initial condition for the evolution of the partonic 
system in ion-ion collisions. The CGC is the source of the thermalization which occurs 
at rather short distance of the order of l/Qs- Our formula is CGC ( saturation) is the 
initial condition for hydrodynamical evolution in A + A collisions.; 

Fig- ED shows that the CGC app roach is able to describe the experimental data using these 
three assumptions (see Ref. |27|). 
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'he CGC approach, based on a new edition of our three assumptions, for deuteron-nucleus 
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pt DISTRIBUTION: PROTON-PROTON COLLISIONS. 



The transverse momentum distribution is a very sensitive check of the form of our input 
for the parton densities. In Ref. 



28[ is shown that Eq. ()13p perfectly describe the pt and ^ 
distribution in proton-proton collision and DIS with the proton target (see Fig. I25p. 
Pt DISTRIBUTION: DEUTERON-NUCLEUS COLLISIONS. 



In Ref. |27[ it is illustrated that the CGC approach is able to describe the transverse 
distributions of the produced hadrons in the deuteron-nucleus collisions n(see Fig. |2ni This 
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FIG. 25: pt and distribution in proton - proton collisions and DIS with the proton target. 



Pictures are taken from Ref. 



28|. 



is a powerful check of our two first assumptions as well as the simplified form of the parton 
densities in Eq. (IT!?j) . 

Pt DISTRIBUTION: ION-ION COLLISIONS. 



In Ref. 



29 



301 1 the good agreement with the experimental data is achieved using the 
CGC theory as the initial condition for evolution which is treated in hydrodynamics (see 
Fig. 1221 and Fig. l28|l . This gives a strong argument in the support of out third assumption. 

AZIMUTHAL CORRELATIONS 

The last issue that I want to touch in this talk is the azimuthal correlations since they 
give a beautiful example of our expectation in the CGC approach. In the new phase-CGC, 
we expect decorrelations mostly because the source of the production of two particles with 
definite azimuthal angle between them is just independent production of two hadrons (see 
Fig. IH^ . Such a production is proportional to the square of parton density and it gives 
a small contribution in pQCD phase of QCD. In this phase the typical process is back-to- 
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FIG. 26: "pt distribution in the deuteron - gold collisions in the CGC approach (see Ref. 



0) 



back correlation due to the hard rescattering of two partons that belong to the same parton 
shower (see Fig. IHUj). 



Fig. |2ni and Fig. ED show our estimates 



12j | of the correlation in one parton shower 



(Fig. and the resulting correlations (Fig. l!?T|l including the production from two parton 
showers. Our estimates certainly reproduce the experimental data of Star collaboration ( 
see Fig. | 



RESUME 



I hope that I convinced you that the distance between experimental data and our micro- 
scopic theory: QCD, is not so long as usually people think. The RHIC experimental data 
support the idea that at high parton density we are dealing with a new phase - Colour Glass 
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FIG. 28: pt distribution in CGC + hydro 
model of Ref. 
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Condensate. The parton density saturation, the most essential property of the CGC phase , 
manifests itself in the RHIC data. The most essential characteristics of the data that can be 
easily and elegantly described in the CGC approach are (i) energy, rapidity and centrality 
dependence of the particle multiplicities; (ii) the suppression of the rapidity distribution in 
the deuteron-nucleus collisions in comparison with the proton-proton collisions in forward 
region of rapidity; (iii) suppression at ?7 = for rather small transverse momenta in the 
same data and (iv) observed azimuthal decorrelation in gold-gold collisions. 

I hope that future experiments at RHIC and LHC will provide more arguments in favour 
of the CGC and our understanding of main features of QCD will become deeper and more 
transparent. 
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